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STABILISATION OF LIQUID METAL ELECTROLYTE SYSTEMS 



Field of the Invention 

The invention relates to liquid m^al electrolyte systems and is especially, though 
not exclusively, applicable to improving the efficiency and reducing the operating 
costs of modem-day aluminium reduction cells. 

25 

Review of the Art currently known to the Applicants 



30 



The invention will be exemplified, and will subsequently be described and 
illustrated in this present specification, with reference to aluminiiun reduction or 
smelting. 
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Modem aluminium production plants consume huge amoimts of electricity. 
Virtually all of them operate by reducing alumina in electrolysis cells or, as they 
are called, pots. In practise a commercial aluminium smelting plant will consist of 
several hundred such pots and will operate on a continuous production basis. 

There are two remarkable features of this process. First, it has remained virtually 
unchanged for over a century since it was first successfully developed (and indeed 
it is still universally known as the Hall-H6roult process after the two scientists 
who first independently discovered it). Second, the amount of energy consumed 
by the process is quite staggering. 

It has been estimated that the modem-day production of aluminium consumes 
about two per cent of all electricity generated woridwide (!) and yet much of this 
energy is absorbed in overcoming resistive losses in the poorly conductive highly 
resistive electrolyte layer of each individual smelting cell. The primary electrical 
driving current can be of low voltage but must be of relatively enormous 
amperage in order for the process to work, given these drawbacks, and it follows 
that any modification which enables that current, the electrolyte thickness, or 
both, to be reduced at all would indeed produce reductions in energy consumption 
which could truly be described as significant in relation to those needed without 
the modification today. 

Efforts have been made, naturally, to overcome this problem but the main limiting 
factor is that, if the electrolyte thickness is reduced beyond a certam critical level, 
instabilities begin to occur at the interface between the liquid electrolyte and the 
liquid aluminium. These instabilities, which manifest themselves as a sloshmg of 
the liquids within the cell, have been the subject of intensive research for some 20 
years or more. In effect, these are interfacial gravity waves, modified by the 
external magnetic fields which pervade the cell and when a certain stability 
threshold is exceeded, these waves can grow by absorbing energy fi-om the 
ambient electric and magnetic fields. 
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The good news is that the wave period is measured in minutes and its growth rate 
in hours, and so the problem ought to be susceptible to some controlled solution. 
The real problem is that once such a wave takes hold, it can disrupt the 
5 electrolysis to such an extent that the cell must be withdrawn from operation. In 

an extreme case, it could destroy the entire cell 



Previously proposed means for trying to eliminate these instabilities include: 



10 • Placing baffles in the aluminium to break up the long-wavelength waves 

whilst relying on friction to dissipate the short-wavelength components. 

• Have a sloping cathode block so that the alimiinium continually drains 
away. 

15 

• Destroy the standing waves by placing hydraulic energy absorbers at the 
edges of the cell. 

• Tilt the anode in harmony with the wave so that the electrolyte layer 
20 remains ahnost uniform and thus one eliminates the perturbation in 

current. 



The first of these prior proposals remains simplistically attractive but both it and 
the second one are limited by the need to find, in practical environments, a 
25 material which survives the chemically aggressive environment in a smelting cell. 

The second option has another difficulty in that thin aluminium layers will not 
properly wet the cathode and this cannot easily or cheaply be overcome. Whilst 
the third option is self-explanatory, the most recent research has concentrated on 
the final one but as far as is known, no practical embodiment has yet emerged. 
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Furthermore, a paper published by Elsevier Science dated 12 November 2001 by 
authors A. Lukyanov, G. El and S. Molokov, is deemed to be relevant to the 
present application as it defines the general background of the mechanism of 
instability, however primarily in the context of determining the reflection 
coefficient rather than proposing a practical solution for controlling instability in a 
cell as is one of the objectives of the present application. 

In summary, despite the length of time the problem has been around, and the 
importance of modem aluminium production to the progress of industrialised 
society as a whole in an era when, paradoxically, conservation of energy is 
becoming more urgent than ever, instability of aluminium reduction cells remains 
the central and unsolved problem in the industry at large. 

The Inventive Concept 

The £q>plicants are proposing a modification of existing current-driven liquid 
metal electrolyte systems (of which an aluminium reduction cell is the obvious 
but not limiting example) which starts fi^om a point quite different from any of 
those outlined above - but which could, we believe, be used in any appropriate 
combination with some, all, or any of the prior proposals outlined above. 

In essence, we impose on such a system an additional, external, magnetic field 
whose design and operating parameters are so chosen as to enable the electrolyte 
thickness to be reduced significantly in relation to those needed without the 
modification. By doing this, we address the very source of the instability, which 
happens due to the interaction of the currents induced by the interface motion with 
the external magnetic field. 

Based on our understanding of the fiindamental mechanism governing the 
instability, we believe it to be possible that, with appropriately designed coils, a 
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ring current around the cell inducing an automating magnetic field will stabilise 
the cell to an appreciable if not total extent. 

Thus, rather than trying to understand fully all the processes happening inside the 
cell we effectively suppress the fluctuations by imposing a suitably powerful and 
time-dependent magnetic field aroimd it. A modem aluminium (or any other 
metal) reduction cell is a complex and highly optimised device. There are a 
multitude of supple physical and chemical processes occurring within such a cell 
and many of them will inevitably interact. A small change in any one parameter 
could well have quite unexpected consequences and these may or may not be 
either inter-related or predictable at all. The size alone of the primary driving 
current makes it ahnost impractical to try to make relatively small adjustments to 
any one aspect of cell operation - for example, the "anode-tilt" approach 
exemplified in the fourth prior proposal outlined previously - with any real 
guarantee of even partial success. 

We by contrast take an overview and we believe that, with appropriate design and 
with tiie ability to adjust the controlling parameters (i.e. field amplitude, 
firequency, and constant background) we are more likely to achieve real 
suppression of instability in a practical format within the foreseeable future. 

In a subsidiary aspect, the magnetic field applied is essentially a vertical magnetic 
field. In this direction, significant effect on instability in liquid metal electrolyte 
occurs which allows the thickness of electrolyte itself to be reduced below levels 
at which conventionally instabihty would occur. 

In a further subsidiary aspect, the magnetic field is dependent on an amplitude and 
fi-equency whose values are approximated through wave reflection analysis on an 
infinite wall. This is advantageous as it allows an appropriate magnetic field to be 
rapidly determined rather than relying on the skilled man to determine an 
adequate field through more extensive analysis. 
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Brief Description of the Drawings accompanying this text 
The accompanying drawings: 

Figure 1 shows diagrammatically an example of a modem Hall-H6roult cell; 

Figure 2 presents the electrolysis zone of the cell schematically; 

Figure 3 shows graphically the existing and the modified instability levels 
occurring in respectively an unmodified and a modified cell in accordance with 
the invention; and 

Figure 4 shows, again in schematic form, one possible set-up embodying the 
invention. 

Figure 5 shows, a schematic diagram of a two-layer system. 

Figure 6 shows, a schematic diagram of wave reflection on an infinite plane wall. 

Figure 7 shows graphically the amplitude of interfacial wave for two electrolyte 
thicknesses when no alternating field is s^plied. 

Figure 8 shows graphically the amplitude of interfacial wave for a reduced 
electrolyte thickness cell with and without the alternating magnetic field. 

Description of the Proposal in Outline 

An example of a modem Hall-Heroult cell generally referenced 1 is shown in 
Figure 1. Cell 1 comprises covers 2, carbon anodes 3, molten salt electrolyte 4, 
molten aluminium 5, collector bars 6, carbon lining 7 and a carbon bus 8. All of 
these components may be of standard kind, modified or substituted if necessary 
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by other relevant components or groups of components by the person skilled in 
the art without any recourse to inventive thought 

The current used in the electrolysis enters the electrolyte zone vertically through 
the anode and is collected by the cathode at the bottom. The thickness of both 
5 layers, electrolyte and aluminiimi, is very small in comparison with the horizontal 

dimensions. Schematically, the electrolysis zone can be presented as shown in 
Figure 2. 

The major part of the consumed energy is wasted in the form of resistive losses in 
the poorly conductive electrolyte, layer 2 in Figure 2. But, when the depth of 
10 electrolyte is reduced below some critical level or the current exceeds some 

critical value, the cell becomes unstable. In other words, the waves at the interface 
between the two liquids start growing. The increment of the resulting instability is 
shown in Figure 3 (curve 1). 

IS It is proposed to apply an external, alternating magnetic field and to regulate the 

currents induced by this field so as to control or even suppress instability. The 
sketch of a possible set-up is shown in Figure 4. In this figure a ring current 
around the cell induces an altemating magnetic field. Jn practice, the altemating 
magnetic field may be created, for example, by coils surrounding the cell or other 

20 means carefiilly selected by the person skilled in the art. The resuh of simulations 

for the circular cell, which exemplifies the most unstable case, is presented by 
curve 2 in Figure 3. One can see that the instability disappears. Analysis of a more 
realistic, rectangular cell shows that the method works successfully in this case as 
well (Figure 8). It is believed that the method may be adapted by the person 

25 skilled in the art for any cell geometry. 

Description of the underlving theorv and exemplar/ results 

In the following description, flow stabilization by an alternating magnetic field is 
30 presented, and the effect of suppression of instabilities is demonstrated on an 

example of a two-layer system in a rectangular geometry. 
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A) Mathematical model of the MHD-modified interfacial gravity wave dynamics 
in a closed domain 

5 Consider the system of two electrically conducting liquids (liquid metal and 

electrolyte) carrying electric current of density J and exposed to a magnetic field 
B presented in Figure 5. 



10 



In the equilibrium state, 

J = Jo = (0,0 -/o), B = (b^B^bJVx[J^ X B J = 0. (1) 



Here (x,y,z) are Cartesian co-ordinates. The last relationship implies that the 
vertical component of the magnetic field Bqz can be arbitrary (given by the 

IS external circuit). 

Let the thickness of the liquid metal layer in the equilibrium state be H\ and that 
of the electrolyte be Any deviation of the interface fi-om the equilibrium state 
(which is inevitably present in a real cell) induces the redistribution of the current 
(and, consequently, of the magnetic field). This process is accompanied by the 

20 wave motion of the two-layer liquid system, hi the absence of the electric current, 

the system is stable (the amplitude of the initial perturbation of the interface does 
not grow in the process of the wave propagation). Eventually, because of the 
natural dissipation in the system, the wave would fade out. In contrast, when the 
cxirrent is on, interaction of the current perturbation with the external magnetic 

25 field can enhance the wave motion and lead to the uncontrolled growth of the 

interfacial wave amplitude. 

The dynamics of the two-layer system is governed by the following equations: 

+ (u, • V)u,] + V(/^. + p,.g2) = Fj - D, 
* (2a) 

30 V.U, = 0, VJ, = 0, V.B = 0, (2b-d) 
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where /= 1, 2 is the layer number of Figure 5; pi is the density; ii, is the fluid 
velocity. Pi is the hydrodynamic pressure, J, is the electric current density in the 
layer (which includes variations induced by the wave motion), B is the total 
magnetic field (which includes the field induced by the external circuit), t is time, 
F/ = J,- X B is the Lorentz force, D, is dissipation describing energy losses in the 
layer. The dissipation term is taken in the conventional form for shallow-water 
equations, i.e. D/ = v,U/, where v/ is the dissipation coefficient. 
The boundary conditions for the two-layer hquid system placed into the poorly 
conducting bath are: 

(u/.ii)^rii=0; (3) 

(JU • ^)side^ls = 0 ; (Ji . n)bottom = -Jo I (Jl . H - J2 . n)i„terface = 0, (4a-c) 

where n is the unit vector normal to a particular surface. 

The current boundary conditions (4) imply the following ranking of 

conductivities: <Tside^b((^2((<^boitom((<^)y which is characteristic for industrial 

aluminium reduction cells (typically, = 33*lO\Om my\ = 200{Om.m)''\ 

=2.10^(Qw.m)-\ a,^,^,^ «0). 

The system of Equations (2), together with the boundary conditions (3), (4) fully 
defines the motion of the two-layer system. 

In the following, the deviation of the interface z = h(x,y,t) from the equilibrium 
state at z = 0 will be discussed. The system of governing equations (2) can be 
significantly simplified if two small parameters are introduced as suggested by the 
actual physical and engineering conditions in the aluminium reduction cells, i.e. 

• 8 = /^i / Z « 1, the shallow water parameter. Here L is the horizontal 
dimension of the cell. Typically, e oc 0.01 . 

• 5 = max hi H\«\, where max h is the amplitude of the interfacial wave. 
That means that we are interested in the dynamics of the small-amplitude 
perturbations, which is perfect for the stability analysis. 
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Implementation of these two parameters means that to the first order in 5 the 
motion of interface is essentially two-dimensional and the following relationships 
are valid: 

5 ni(x,y,zj) «8v/ (x^yj), h {x,y,z,t) »5ii {x,y,t) , Fj {x, y. z, t) « Sf, (x^yj), (5) 

where V/, t^, f/ are new, unknown, (9(1) functions. These are the normalised 
velocity, and the perturbations of the interface and the Lorentz force, respectively. 

10 Taking into account the shallow-water, small amplitude approxunation (5), the 

analysis of the original equations (2H4) shows that to the first order in <5 the 
following conclusions can be made: 

• the current perturbation induced by the interface motion is horizontal, i.e. 
J« Jo + j||(^*>.0 (here and elsewhere subscript || denotes a component of a 

1 5 vector in the (x,3;)-plane), 

• the Lorentz force acting on the liquid metal depends only on the vertical 
component of the external magnetic field: fi » Boz, 

• The Lorentz force acting on electrolyte is much less than that on the liquid 
metal. i.e. |f2|«|f,|. 

20 

As a result, one can conclude that by controlling the vertical component of the 
magnetic field Boz (which is givra by the external circuit) it may be possible to 
control the force inducing the unstable motion of the interface. One such a 
possibility is to superpose a certain alternating magnetic field onto the external, 
25 stationary field. 

Therefore, the following form of the vertical component of the field is considered: 

where Bqz is a constant, while fimction b(x, y, t) can be arbitrary. In previous 
30 studies, the magnetic field has been supposed to be stationary (i.e. independent on 

time) and fixed. 
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Under all assumptions made above the system governing motion of the interface 
assumes the form 

0-c*V^T, = c^V<|).[Vx6(x,3'./)cJ-v,$. (6) 

ot dt 

^'<l>=-Pri, (7) 

Here V.e,%^ + e,^; c^4(P^-9^)^J H.-^Pz' H^V is the speed of the 
inter£acial gravity waves in the absoice of the external magnetic field, 
(|)(x.;;,0 = a,5„g-'(p, -p,)-'q)(*,:v.O normalised electric potential (i.e. 

h = -<'^^<P ), and P = -^"^o 4^.^2(P, - ?2)g] 



It should be noted that natural dissipation in the cells plays crucial role in the 
stability of the existent set-iq)s. Typical value of the non-dimensional parameter p 
in this case is ^20. Without dissipation the stable operation is only possible for 
small values of pa 1 which are totally impractical. 
1 5 The boundary conditions (3), (4) yield: 

^ = 0. ^ = -bix,y,0^ atr (8a,b) 
on on OT 

Here the function r(x, y) = 0 defines the shape of the boundary (horizontal 
geometry of the cell); dy^^and stand for normal and tangential derivatives 
tor = 0, respectively. 

20 Analysis of the system of Equations (6) - (8) by the skilled man in the art in the 

simplest case, when b ^1 (uniform, constant magnetic field), has revealed the 
mechanism underlying the interfacial mstability. Essentially, it has been shown 
that the instability (if it occurs) is inspired at the boundaries of the cell by the 
wave reflection with the reflection coefficient greater than 1. Earlier studies 

25 missed this very point of the instability mechanism for a uniform external 

magnetic field. For this type of fields the first term in the right-hand side of 
Equation (6) vanishes, and Equation (6) becomes essentially decoupled fi'om 
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Equation (7). It is the boundary condition (8b) that is responsible for the 
development of the mstability. And here is the remedy: there is an arbitrary 
function b(x, y, t), which is essentially the externally applied magnetic field, in 
this boimdary condition. 

Derivation from this underlying theory enables the preferred inventive external 
magnetic field b(x, y, t) to be found which leads to the attenuation or even 
suppression of instability. 

Below results are presented for the simplest case of spatially uniform alternating 
magnetic field 

i = l+Z>oCos(fi)or + ^o) (9) 

Here bo is the normalised amplitude, wq is the fi-equency, and 6o is the initial 
phase of the controlling external magnetic field which is to be obtained. 

For a realistic geometry of the cell the problem defined by Equations (6)-(8) must 
be solved numerically. For calculations in the specific case of a rectangular cell as 
presented hereafter, second order central differences may be used throughout. 
Equation (6) may be discredited using an exphcit scheme in time. A fast Poisson 
solver may be used to solve Equation (7). 

For calculations 32 points per unit length may be used. The scheme has been 
successfully tested on several benchmark problems to ensure both high accuracy 
and the absence of numerical dispersion. Other methods of determining 
advantageous magnetic field types may be employed and will be selected by the 
person skilled in the art firom known alternatives. 

An approximation for parameters bo and (Do can be advantageously 
obtained from the corresponding problem of a reflection from an infinite plane 
wall (see Sec. B). Starting from these initial estimates the frequency and the 
amplitude are either increased or decreased to minimize the increment of 
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instability. The parameters are adjusted iteratively until stability of the interface is 
achieved. 

B) Approximate determination of the amplitude and frequency of the external 
magnetic field: reflection from the infinite plane wall 

One example of reflection analysis from an infinite plane wall is presented in this 
section. 

Both the amplitude and frequaicy of the controlling parameters of the external 
magnetic field are estimated using the simplest model of the reflection of the 
plane wave from the infinite boundary m the absence of dissipation as shown in 
Figure 6. 

In a previous study of this kind where b was assumed si, the reflection 
coefficient \i was found to be greater than 1 for some angles of incidence. In other 
words, the wave was being amplified at the boundary. It is clear that in the 
presence of the alternating magnetic field b(t) given by Equation (9) one obtains 
|i=ji(6fttt)ft;. Now we are going to find such controlling parameters bo and oo that 
the reflection coefiBcient n^l. To achieve this it is convenient to represent the 
problon of the reflection of the plane wave from the wall in the form of the 
integral equation for they-Fourier component of jj (x, y, t). 

The dependent variables in the reflection problem are represented in the form 



where ky is the wave number of the incident wave. 

The Fourier-transform with respect to x leads to the following integral equation 
for the function r}{x,t) at the boundary: 




.3^ 
' dx 




(10) 
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while function ^(^»^) satisfies equation 

-c^^ + A:VJ7=0. (11) 



Applying further the Fourier transfonn over t i.e.7(jc,r)= J77^exp(-ifi)^)rfa>, 
gives the solution of Eq. (1 1) as follows: 

'7e> = Q(^»)exp(i*,:r) + C2(fi>)expHV)^ (i2) 

where =^o)^fc^ -k^ , and Cy{a)\C2(o)) are spectral powers of incident and 

reflected waves. Substituting Equation (12) into Equation (10) yields a functional 
equation, which links spectral powers of the reflected wave and the incident one, 
namely 

' 0,(0)) , C,(Q)) 1 



*,(o)){C,(a))-Q(o>)}. 



C,((0±0),) ^ C,(CQ±Q),) ] ^ 

°[^^+/^,(©±0)o) i^-z^,(co±a)jj 



(13) 

Equation (13) can be solved iteratively assuming the spectral power of the 
incident wave is given, for instance Cj(a)) = 1. This gives the values of 6^ and o)©, 
which can be used as a starting point in our analysis of the instability in the 
rectangular cell with dissipation. Further, these parameters must be tuned using 

15 developed numerical code to achieve stability. 

It is worth noting that Equation (10) can be used to solve a more general, inverse 
problem. That is, if one prescribes the spectral power of both the incident and the 
reflected waves, then one can obtain necessary time dependence of the controlling 
magnetic field b(t) rather than assuming any parametric form of the type (9) a 

20 priori. 



C) Control of ingtahili tv in a rectangular cell 



1* 
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The stabilizing effect of an alternating magnetic field in a rectangular cell will be 
demonstrated on the following example. Let the geometrical parameters of the 
cell be: the length L\= 9.8m, the width L2- 3.4m, the thickness of electrolyte layer 
H^-Scm and the thickness of aluminium layer H\- 25 cm. The total current 
5 passing through the cell is /c = 175 kA. The constant external magnetic field has 

been taken to be Bo = 3J(r^T. These conditions correspond to a stable process of 
the aluminium production, which is confirmed by computer simulations and 
corresponds to horizontal curve in Figure 7. 

If one reduces the thickness of electrolyte layer by 5 %, i.e. Hz^JS, the cell 
1 0 becomes very unstable. Such resulting instability is shown in the growing curve of 

Figure 7. As is seen, the growth rate is rather significant and after 30 minutes 
short-circuits occur. 



Figure 8 shows a stabilized cell with reduced thickness of an aliuninium layer by 
1 5 an alternating field. 

Operation of the cell with the same (reduced) electrolyte layer thickness but with 
application of an alternating magnetic field is given by Equation (9) where bo = 
0.66, coo = 20 radsec'^ 9o = 0 is shown in Figure 8. 

The proper fi-equency ©o and the amplitude 60 have been found according to our 
20 method described in section B) and tuned to achieve stability. Starting values 

were not far away from those that give stable operation, i.e. 
^cpprax (ol'P"^ ^AQrad sec-\ Note that bo is normalised with Bo = 

As a result the cell becomes stable as Figure 8 shows. As will be shown in the 
25 following section, this result is particularly encouraging in terms of actual energy 

savings. 



Reduction of the energv consimiption 



30 



Let's calculate energy losses in the electrolyte layer per one millimetre under the 
parameters listed above. The conductivity of molten electrolyte is Gc == 200 



wo 03/057945 



16 



PCT/GB03/00072 



(Om.m)'\ Then in each millimetre of the electrolyte layer (AL = 1mm) energy 
losses due to Joule dissipation areiFTg = llAL/{ag L^L2} = ^.6kWatt , Since the 
inventive magnetic field application has permitted the electrolyte layer's thickness 
of being reduced by AHj = 2,5mm , it follows that the electric energy 
consumption may be reduced by AFF^ = 1 1 .SkWatt . On the other hand, to create 
the stabilising external, alternating magnetic field by a coil, one needs to spend no 
more than = 57Wati^ provided the coil has 300 loops of copper wire of 0.5 cm 
in diameter. 

So, the ratio is just WjAW^ =0.5%. That is, the energy expenses for the 

production of the controlling magnetic field are very small in comparison with the 
resulting savings. 

Two-layer systems carrying electric current in the presence of a magnetic 
field can be stabilized by the £q>plication of an ext^nal altemating magnetic field. 
Calculations for a typical geometry of an industrial aluminiimi reduction cell in 
the presence of a uniform field show that the energy losses required for 
stabilization are minimal. 

Similar calculations may be performed for cells of various shapes and even for 
spatially non-unifoim magnetic fields. The person skilled in the art will adapt the 
preceding broad underlying theory for each specific case while the current scope 
of the invention is defined in the claims which now follow. 



